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Anomalous Relations between H and Z Components of 
Transient Geomagnetic Variations* 


By Tsuneji RIKITAKE and Izumi YOKOYAMA 


(Earthquake Research Institute, Tokyo University) 


Abstract 


This is to summarize the writers’ studies on the anomalous 

behavior of geomagnetic variations in Japan which were published in 

a series of papers in the Bulletin of the Earthquake Research Institute 

[1]. An anomalously large amplitude of the vertical component of 

geomagnetic variations is pointed out at Hermanus, S.A., and some 

observatories in Japan. The distribution especially in Japan is statis- 

tically examined. The anomaly is found in a limited area in the 

central part of Japan. The results of the analysis based on the 

potential theory show that the anomaly is caused by the magnetic 

field originating in the earth. The anomaly is well approximated by 

‘an apparent magnetic dipole situated under central Japan. But it 

still remains unknown why such a localized anomaly occurs on occa- 

sions of geomagnetic variation. 
- 1. Introduction. 

The vertical component (4Z) of geomagnetic variations of short period is 
generally small compared to the horizontal component (4H) as we can see on 


magnetograms from various observatories. In some places, however, we find an 


aZ 4Z 


Fig. 1 The relation between 42 and 4H for Fig. 2. The relation between AZ and AH for 
short period variations during the short period variations during the 
magnetic storm on Aug. 3, 1949 at magnetic storm on Aug. 3, 1949 at 
Hermanus. Kakioka. 


* Dedicated to Professor M. Hasegawa on his sixtieth birthday. The writers are greatly obliged 
to him for his ingenious method of graphical integration which is used in this paper. 
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anomalously large 4Z. For instance 4Z/4H amounts to as much as 0.83 on an 
average at Hermanus Magnetic Observatory, South Africa, a remarkable parallelism 
between H and Z-curves being also observed there. The relation between AZ and 4H 
at Hermanus is illustrated in Fig. 1 with respect to the variations during a magnetic 
storm that occurred on Aug. 3, 1949. Meanwhile 42/4H amounts to 0.56 on an 
average at Kakioka Magnetic Observatory, Japan as also shown in Fig. 2 for the same 
storm. According to the writers’ detailed investigations [1] concerning the distribution 
of 4Z/4H in Japan, the fact that the anomalous behavior of geomagnetic variations 
is confined to the central part of Japan has been clarified. The shape of Z-curves 
observed at a station along the coast of the Japan Sea, only 200km distant from 
Kakioka, is quite different from those observed at Kakioka. Even the sign of 4Z 
reverses at stations in the northern part of Japan, while 4H seems to be fairly regular 
throughout the country. The purpose of this paper is to study the anomalous 
behavior of geomagnetic variations and its possible cause by summarizing the various 


studies in the writers’ previous papers [1]. 


2. The distribution of 
4Z/AH in Japan. 

In order to bring 
out more clearly the 
anomalous’ characteris- 
tics of geomagnetic 
variations in Japan, the 
writers collected mag- 
netograms from perma- 
nent and temporary 
observatories in Japan 
and its vicinity. Fortu- 
nately we had well-dis- 
tributed stations there, 
especially on the occa- 
sion of the magnetic 
observation related to 
the solar eclipse on June 
19, 1936. In addition to 
those observations the 
writers added two sta- 


Fig. 3. The distribution of the statistical value of AZ /AH in tions for continuous 


Japan. The variations in H and Z com i 

ponents at recording of three - 
the time of the bay at 23h on Sept. 11, 1952 are also = ik 
illustrated for respective stations. magnetic Compe oe 


Aburatsubo (139°37'E, 
Although the periods of these observations 
the characteristics of geomagnetic variation seem to be nearly 


35°09’N) and Maze (138°48’E, 37°44/N). 
are not the same, 
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invariable. Hence we may construct a chart in which the distribution of the statis- 
tical values of 4Z/4H are illustrated as shown in Fig. 3. In the figure, the records 
of a geomagnetic bay that occurred at 23h on Sept. 11, 1952 are also shown for 
respective stations. We can clearly see that the amplitude of 4Z is very large near 
the central part of Japan. The writers have never heard of such a strange distribu- 
tion of geomagnetic variation except for the opposition of sign of 4Z in British Isles 
and Western Europe [2]. 

3. Anomalous geomagnetic fields originating within the earth. 

As the first step in clarifying the origin of the strange distribution mentioned 
above the writers intended to separate geomagnetic changes for each station into the 
parts originating outside and inside the earth. Since we had many; temporary stations 
in and around Japan at the time of the solar eclipse on June 19, 1936, the writers 
collected copies of magnetograms of the magnetic storm that occurred with a SC at 
9h 41m on June 18 from various observatories all over the earth. On the basis of 
these data the writers calculated the distribution of magnetic potential W for the 
SC and a remarkable change that occurred at 5h 50m GMT on June 19. The method 
of the graphical integration which was devised by M. Hasegawa and M. Ota [3] was 
used in the calculation. The equi-potential lines are respectively illustrated in Figs. 4 
and 5 together with the equivalent overhead current arrows. Combining the distribu- 


japan 


Bi) 
sz = 
ae >) Wk * ia nner 


Z 


Bi 
14s 
oe 

on) 
aes 
| tss 
\ 
au : 
MS 
ES 


IN a7 


is 


are 


Py he 


503 


Fig. 4 The distribution of magnetic potential for the SC at 9k 41m GMT on June 18, 
1936 as calculated by means of the graphical integration method. T he unit of 
the potential should be read in gamma multiplied by the earth’s radius. The 
equivalent overhead current arrows are also shown for respective stations. 
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Fig. 5 The distribution of magnetic potential for the variation at 5 50m GMT on June 


19, 1936 as calculated by means of the graphical integration method. 


The unit 


of the potential should be read in gamma multiplied by the earth’s radius. The 
equivalent overhead current arrows are also shown for respective stations. 
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Fig. 6 The equi-potential lines for the 
external origin part of the magnetic 
potential for the SC. The unit is 


gamma multiplied by the earth’s 
radius. 


180° 195° 


The equi-potential lines for the 


internal origin part of the magnetic 


potential for the SC. The unit is 


gamma multiplied by the earth’s 
Tadius. 
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Fig. 8 The distribution cf the external part 
of the magnetic potential for the 
variation at 54 50m GMT on June 


19, 1936. The unit is gamma 


multiplied by the earth’s radius. 


210° 

Fig. 10. The distribution of the external part 
of AZ for the same variation. The 
unit is gamma, 


The distribution cf the internal part 
of the magnetic potential for the 
same variation. The unit is gamma 
multiplied by the earth’s radius. 


180° 195 210° 
Fig. 11 The distribution of the internal part 
of AZ for the same variation. The 
unit is gamma. 


tion of W with that of 4Z, the external and internal parts of W and 4Z are calculated 
for stations in and around Japan by means of E.H. Vestine’s method of surface 


integration [4]. 


Owing to the irregular distributions of 42, the accuracy of the 


calculation is not too reliable in general. But as far as the stations ia Japan and its 
neighborhood are concerned, the results should be approximately right because we had 
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fairly many stations there. Both parts of the potential for the SC are shown 
respectively in Figs. 6 and 7, while similar results are illustrated in Figs. 8 and 9 for 
the change at 5/2 50m on June 19. The external and internal parts of 4Z for the 
latter are also shown in Figs. 10 and 11. As can be clearly seen in these figures, the 
distribution of the external origin parts are fairly regular for both variations. On 
the contrary, the internal origin parts are distributed in a strange way. 

Although the anomalous distribution of the internal origin part can be seen less 
clearly in the case of the SC partly because of the smallness of the amplitude of the 
variation, we find some characteristic equi-potential lines as well as equal variation 
lines for 4Z of internal origin concentrated near the central part of Japan. If we 
consider that the variations of internal origin are caused by electric currents flowing 
just under the earth’s surface, the currents should flow in the directions which are 
shown in the figures. It seems quite strange that the induced currents flow to the 
same direction as that of the inducing currents provided the usual standpoiat of 
electromagnetic induction within a homogeneous earth is adopted. Hence we must 
imagine some particular conditions beneath Japan. 

In Figs. 10 and 11, we can see that the external parts of 4Z direct upward 
while those of internal parts direct downward taking almost the same magnitude in 
the northern part of Japan. This seems to be the reason why the amplitude of 4Z 
is usually small there as stated in the last section. But 4Z of the internal origin is 
definitely large and positive near the central part of Japan, this being the reason 
why it overcomes the external part causing a large 4Z of the positive sign there. 
Taking into account these tendencies of external and internal parts of geomagnetic 
variation, the strange distribution of 4Z/4H as shown in Fig. 3 can be phenomeno- 
logically explained. 

4. The origin of the anomalous behavior of geomagnetic variations in Japan. 

As is made clear in the last section, the anomalous behavior of geomagnetic 
variations in Japan is caused by magnetic fields originating in the earth. However 
the internal origin part of geomagnetic variation contains the part produced by 
electric currents induced in the electrically conducting region the depth of whose upper 
surface is estimated at about 400km [5, 6, 7, 8}. In order to abstract the anomaly 


itself, the part originated in this region should be eliminated. On the basis of one 
of the writers’ studies (T.R.) [8], this elimination was carried out for the variation at 
5h 50m on June 19, 1936. After the elimination we find a particular distribution of 
the magnetic field in Japan. This field is well approximated with that of a magnetic 
dipole situated at a depth of 150km under the central part of Japan, the direction of 
it being approximately north and having a dip of 20°. Though this expression with 
a dipole is nothing but a conventional model, we must pay attention to the fact that 
the anomaly is well explained by considering a localized source. We can say nothing, 


however, why the apparent magnetic dipole occurs on occasions of geomagnetic 
variations, 
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5. Conclusions. 

The writers pointed out an abnormally large 4Z of Short-period geomagnetic 
variations observed at Hermanus as well as some Japanese observatories The 
characteristics in Japan were statistically examined with data from well -distriboted 
observatories. Meanwhile it was made clear that the anomalous behavior in Jap=n 
is caused by the magnetic field originated within the earth by anelyring world wide 
data from the standpoint of the potential theory. Detailed reports were already 
published in a series of papers [1] in the Bulletin of the Earthquake Research Instimte 

In conclusion, the writers are grateful to Professor M. Hasegawa, Dr. M. Ora, 
Professor T. Nagata and Mr. N. Fukushima who are interested in this study and offer 
useful data to the writers. The writers are also indebted to many geomagneticizns 
in foreign countries who kindly sent copies of megnetograms at the writers” request 
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The Thermal Fluctuation After Effect found in the Natural 
Remanent Magnetic Polarization of Rocks 


By Naoto KAWAI and Shdichi KUME 


(Geological and Mineralogical Institute, Faculty of Science, Kyoto University) 


Abstract 


In the study of natural remanent magnetic polarization of 
rocks, the authors found many kinds of sedimentary rocks whose 
major part of the polarity may probably be attributed to thermal 
fructuation after effect resulted under the influence of the geomagnetic 
field during the geological interval. In this paper are reported the 
- polarity of tertiary rocks observed, and experiments of the related 
problems carried out in our laboratory. 
I. Introduction. 
A few years ago, J.W. Graham” proposed many ingenious methods for examin- 
ing the stability of natural remanent magnetic polarization of rock (abbr. N.R.M.). 
From his numerous field observations it was well confirmed that many rocks kept 
their initial direction of polarization unchanged through the geological time up to the 
present. In order to test the said stability of rocks in Japan, the present authors 
employed one of the Graham’s methods. Many pieces of pebbles were sampled from 
a conglomerate of a known geological age, and their directions were measured by a 
sensitive astatic magnetometer.2> Favorable agreements with the Graham’s result 
were obtained in the case of many igneous pebbles. However, the data from sedi- 
mentary pebbles obviously showed themselves as the exceptions. Those conglomerates 
whose pebbles had their direction of the N.R.M. nearly parallel to that of the geo- 
magnetic field have been found from the tertiary basins in Bésd and Osaka. It is 
quite unlikely that the pebbles in the conglomerate possess such a particular 
direction in the field, when we assume that the N.R°M. is stable through the geological 
time. This seemed to suggest the following presumption that the direction of the 
pebbles might have changed into the present configuration of the polarity from the 
initial and probably random distribution with which the pebbles had been migrated 
into the stratum. 
After the field examinations, each pébble specimen has carefully been kept. in 
the laboratory in constant orientation since these three years and measuring of the 
polarization was successively repeated with proper duration of intervals. A distinct 


change which is exceedingly larger than measuring errors, was observed in both the 
direction and the intensity of the polarity, 
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In this paper the data are summarized with an aim of making connection to the 
Néel’s theories**© on the thermal fluctuation after effect. 


II. Field observations. 


The examinations of the stability were carried out at thirteen localities, five 
sheets of conglomerates being selected from Osaka, as well as eight sheets from B6dsd 
sedimentary basin. About 20 pieces of pebbles were sampled from each conglomerate 
at each locality and their directions of N.R.M. were measured by the astatic magneto- 
meter. The results are illustrated in Fig. 1A-Fig. 1E. with respect to the following 
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A: Igneous Pebbles B: Sediments C: Sediments D: Sediments E: Sediments 
xX: Position of present ©: Positive magnetic pole @: Positive magnetic pole 
geomagnetic field on upper hemi-sphere on lower hemi-sphere 
Figs il 


typical cases. In these figures the positive magnetic poles of pebbles were adopted 
and plotted in the Schmidt’s equal-areal projection. The solid and hollow circles in 
the figures represent the poles which are situated respectively on the lower and upper 
hemi-sphere of the projection and the cross marks correspond to the position of the 
present geomagnetic field. 

The pole distribution of igneous pebbles in a”Miocene conglomerate in Osaka 
group is shown in Fig. 1A. (The age of the bed is estimated from palaeontrogical 
criteria to be older than 10’ years, whereas the mean intensity of the N.R.M. is 
2x10-°C.G.S./cm*). The poles are distributed in random manner on both the upper 
and the lower hemi-sphere of the projection with nearly equal population. The result 
agrees well with the Graham’s observation. Meanwhile, the similar examinations were 
carried out with respect to the sedimentary pebbles in Miocene conglomerate in Boso 
group (The age is nearly of the same order as that of Osaka group, but the intensity 
of the N.R.M., is weaker than the former, ranging from 10-* to 10-°C.G.S./cm‘). 
Throughout the results of the sedimentary pebbles (Fig. 1B—Fig. 1E), the poles of the 
polarization are distributed only on the lower hemi-sphere of the projection. As will 
be seen in Fig. 1D and Fig. 1E, the direction of each pebble is approximately parallel 
to the present geomagnetic field vector. In contrast to the diverging distribution of 
the poles on the projection of Fig. 1A, the directions of these pebbles remarkably 
converge to the cross: mark on the diagrams. 

III. The change of polarization in the laboratory. 

The pebble specimens which had been prepared for the stability test were 
preserved in the laboratory since these three years. Each of the pebble specimens 
was carefully kept with the fixed orientation so that the earth’s magnetic field might 
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apply in the direction appropriately deviated from that of the polarity. And the 
direction was measured successively, allowing the interval duration of 2, 16, 43, 110, 
380, 785 and 1000 days respectively. In almost every sedimentary specimen except in 
the case of the pebble A, changes of the polarization were sufficiently detectable by 
our measuring equipment. The angular change of the polarization is illustrated by 
the pole movement in Fig. 2-3, where the durations of time (in days) are shown by 
the number suffixed to the poles. The following tendencies may be to be noticed. 


\ 
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S 
Bigsene Bigs 

xX: Position of present ©: Positive magnetic pole @: Positive magnetic pole 

geomagnetic field on upper hemi-sphere on lower hemi-sphere 


(1) During the change the direction of the polarization is generally restricted 
on a definite plane which is determined by the two direction, one being the same as 
the initial direction of the polarity before the change, and other that of the geo- 
magnetic field in laboratory. 

(2) The velocity of the pole movement was greater at the beginning of the 
change,” and then has decreased gradually to the present. Within our experimental 
time interval ranging from 2 to 1000 days so far, the angular changes are approximately 
proportional to log ¢. Amount of the angular change vs. log ¢ curves are shown in 
Figs. 4A-4E. The letters A, B, C, D, E in these figures correspond to the curves 
obtained from the pebbles A, B, C, D and E respectively. 
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*— Deviation from the Direction of the Present Geomagnetic Field — 
Fig. 4 
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(3) Although the slopes of the curves take various values with respect to rock 
facies, there may be a qualitative but close correlation: the pebbles with better 
convergence in the conglomerate (in situ) have faster rate of change in the laboratory. 


Comparison of Fig. 1 to Fig. 4 demonstrates this correlation. 
IV. Temperature dependence of the change of polarization. 


Another trial to observe the change of the polarization was carefully examined, 
keeping the temperature of the specimens appropriately constant throughout the 
experimental time interval. One of the pebble specimens was kept in a constant 
orientation so that the geomagnetic field may be applied in the direction perpendicular 
to that of the specimen polarity. And the polarity was measured by the astatic 
magnetometer. The change of the polarity in both its direction and intensity is 
illustrated by the vectorial representa- 
tion of the polarity during the experi- 
ment (Fig. 5). Meanwhile, the pebble 
specimens were ground into powder or 


69.3 x 10-°C.G.S. 


were freely released into grain with 
least mechanical shock as possible. 
The powder was sealed in a 
glass cylinder so that each grain may 
cancel their polarization with other, 
and the test specimens with non- 
permanent polarization were prepared. 
© The specimens were kept during 200 
hours at a constant temperature in a 
thermostat under the application of a 


Applied Field constant magnetic field. By these 


onion treatments the sufficiently measurable 
polarizations were induced in the specimens. And by measuring the intensity with 
proper duration of intervals it was found that the intensity increased with time in 
a manner as is illustrated in Fig. 6. Similar examinations were carried out with 
various constant temperatures, 28°C, 40°C, 55°C and 68°C. The results obtained at 


each experiment are shown respectively © 
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in Fig. 7. 2 

From the experiments it may be ° 
noticed that the intensity of polariza- 
tion induced by the experimental 
procedure is nearly of the same order 
as that of the N.R.M. which we could 


originally find in the field, although the 


5 10 x10 *sec, 
ratios of the former to the latter take ping 


various values in different specimen Fig. 6 


70 N. Kawai and S. Kume 
3 pieces. 

"2 The direct correlation of our 
4 


results to the Néel’s theories on the 
thermal fluctuation after effect** and 
the related experiments by many 
authorities®***® seems to be difficult in 
the present stage. The problem will 
be postponed until the future when 
the sufficient data of our experiments 


Lo 


will have been accumulated. 
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On Distribution of Nitrogen in the Upper Atmosphere* 


By Teruo SATO 
(Geophysical Institute, Kyoto University) 


Abstract 


In this paper the dissociation of the nitrogen in the upper 
atmosphere is studied. It is considered that the nitrogen atom is produced 
by dissociative recombination of the positive nitrogen ion and the 
electron, and by the pre-dissociation by the absorption in the Lyman- 
Birge-Hopfield band. The former is effective in the region above about 
120km, while the latter predominent below that level. The disappear- 
ance of the nitrogen atom’ is mainly due to the recombination by two 
or three body collision. Approximately, the fomer mechanism of . 
disappearance seems to be preponderant over the latter in the region 
where the dissociative recombination is conspicuous. Hence, two 
pairs of equations in an equilibrium state are solved. In the calcu- 
lations, it is assumed that the temperature is constant throughout 
the region; the ratio of the concentation of the oxygen atom to that 
of the sum of the nitrogen is one-half; and the distribution of total 
nitrogen follows the law of gas in a static equilibrium, though the 
distribution of nitrogen molecule with height is not given previously. 
The results show that the nitrogen begins to dissociate at the level 
from 120km to 130km, and proceeds to completion at about 160km- 
220km. The concentration of the nitrogen atom has peak at about 
140km-150km altitude and its magnitude is of the order of 10'°/c.c.. 
$ 1. Introduction 
Although the dissociation of the oxygen has been studied by many workers (for 
example, Rakshit [1], Penndorf [2] and Moses and Ta-You Wu [3)), the dissociation of 
the nitrogen has not been researched in detail. This problem is very important because 
the formation theory of the ionized region, especially, of that above the FE region has 
a close relation to the kind and distribution of the absorbing particle of the solar 
radiation. Since the measurement of the pressure and that of the density in the upper 
atmosphere by the rocket flight have been made, it becomes necessary to know dis- 
sociations of the oxygen and the nitrogen for deducation of the concentration of the 
particle and the temperature in any height. The Rocket Panel [4], for example, 
assumes that the oxygen and the nitrogen begin to dissociate at the levels of 80km 
and 120km, and proceed to completion at 120km and 220km respectively. There are 


* Contribution to Geophysical Papers dedicated to Prof. M. Hasegawa on his sixtieth birthday. 
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two main mechanisms of the dissociation of the nitrogen. One is the dissociation 
suggested by Bates and Massey [5], Bates [6], and Mitra [7], in which the nitrogen 
molecule is ionized by radiation with a wave length below 795 A, the positive nitrogen 
ion produced combines with an electron, and it dissociates into two nitrogen atoms. 
The other is the pre-dissociation suggested by G. Herzberg and L. Herzberg [8], in 
which the nitrogen molecule is excited by radiation with the wave length 1150 ae 
1250A and later this excited molecule dissociates into two atoms. The former 
suggestion is useful in the upper inoized region. Recentely Bates [6] and Deb [9] 
calculated the distribution of the nitrogen atom using these two mechanisms. Bates 
computes mainly the distribution of nitrogen atom at the level of the £ region, 
including the disappearance of the nitrongen atom due to the existence of the oxygen 
atom and the nitric oxide. According to Deb, the dissociation of nitrogen does not 
exceed fifty percent even at the level of 400km. In his study, however, the distribu- 
tion of the nitrogen molecule with height in an equilibrum state is assumed, which 
in computation must be treated as unknown. Hence it appears that the results are 
unsatisfactory. In our study the two mechanisms of dissociation are still adopted, 
and for the disappearance of the atom the recombination by two or three body 
collisions are used corresponding to the two dissociation mechanisms, as mentioned 
in Appendix I. The distribution of the nitrogen molecule with height is not given, 
but it is assumed that the sum of the total nitrogen particles distributes in such 
a way as a gas in Static equilibrium. The ratio of the concentration of the total 
nitrogen to that of the oxygen atom is taken to be 2. 


§ 2. Production and Disappearance of Nitrogen Atom 


(i) Production of Nitrogen Atom 
Bates and Massey, and Mitra suggested that in the F region the following 


ionization and dissociative recombination are effective, 
N.+hv- N.* +e, (a) 

1 
N2*+e>2N. (b) Se 


The energy which ionizes the nitrogen molecule is that above 15.6 eV (795 A). On the 
other hand, G. Herzberg and L. Herzberg suggested that the nitrogen is also dissociated 


by pre-dissociation by the radiation in a range 1150 A-1250A as follows, 
NAX'E,*)+hy > Ni(a'II,), (a) ) 
N,(a'II,) > 2N. (b) * 


According to Bates, this process is effective in the lower region than the F region. 
Hence in the computation, (1) or (2) may be neglected corresponding to the height 
concerned. Bates also pointed out that the nitrogen atom is produced by dissociation 
of NO as follows, 


NO(X?1I)+hv(6.5eV) + NO(A2S* or Cs*), | 


NO(A2S* or C23*) > N+O, (3) 
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This appears to be effective in ‘the same region as™that in which the process (2) is 
preponderant. According to Bates, when (3) and (6), (7) in § 2 (ii) are taken into 
account, the concentration of the nitrogen atom becomes smaller than that in the 
pure nitrogen atmosphere. But even for the latter, the concentration is very small 
compared with the total particle. As we are interested mainly in the upper region 
above the E region, and the computation is complicated we do not take into account 
these reactions. Hence, the results in the lower atmosphere will show the upper limit 
of the concentration of the nitrogen atom. According to Appendix I (a), the range 
in which the process (1) is usable, is above the region where the concentration of the 
total particle is at least 10'/c.c., and process (2) is usable below that level. 

(ii) Disappearance of Nitrogen Atom 

Moses and Ta-You Wu and Deb showed that oxygen and nitrogen atom dis- 
appears by recombination by two-body collision rather than three-body collision, except 
in the very low atmosphere. This reaction for nitrogen atom is given by 


N+N —> Ne+hv. (4) 


But Bates, Penndorf and Mitra showed that the rate constant by three-body collision 
is larger than that by two-body collision, and is of the order of 10-%cm‘*/s. Two 
ways of three-body collision are given by 


N+N+M > N2+M, (5) 
N+0+M > NO+M, (6) 


where M is a third body which takes away the excess of the energy and momentum. 
Another process is given by Bates as follows, 


N+NO —> N,+0. 4 CF) 


We neglect (6) and (7) for the reason mentioned above in (i). In Appendix I (b) it is 
shown that the process (4) is predominent above the level which has the con- 
centration of the total particle of at least 10''/c.c.. Thus the process (1) and (4) 
are considered to be active in the same range of height and (2) and (5) in another 
range of height. 


§ 8. Absorption of Solar Radiation 


(i) The solar energy in a spectral range below 795 A is not only absorbed by 
the nitrogen molecule, but also absorbed by the oxygen atom and the nitrogen atom. 
The oxygen atom is ionized by the radiation with the energy above 13.6 eV (910 A) 
and the nitrogen atom by the radiation above 14.5eV (850 A). According to Bates 
and Seaton [10], the absorption coefficient is 2.6x10-'%cm? for the oxygen atom and 
9x10-'%cm? for the nitrogen atom. Nicolet [11] pointed out that the absorption coef- 
ficient for the oxygen atom is 4.5 10-'%cm? for the radiation between 13.6 eV and 17 eV, 
1.1x10-'%cm? between 17 eV and 18.6eV and 1.6 x 10-'"cm? above 18.6eV, and 2x 10-1!%cm? 
for the nitrogen atom. These values are comparable with the absorption coefficient of 
the nitrogen molecule given by Bates and Massey which is about 1x10-!%cm?. Thus the 
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radiation decreases more rapidly in the case of the absorption by sum of these 
different kinds of particles than the case of the absorption only by nitrogen molecules. 
The concentration of ionized nitrogen molecule produced by the process (1) is 


given by 
q. = Aw,S,0exp[— j (AvLNo]+ AyLN 1+ AolO}}sec x dhJLNo], (8) 
ah 


where brackets represent the concentration of the particle per c.c., A, the absorption 
coefficient referred to the particle 7, which is considered to be independent of the wave- 
length of the radiation, S,.. the number of mean total quanta available, # the coordinate 
of height and x the zenith angle. If we take Ay.~Ay»~Ay=107"cm?, and assuming 


[O]=1/2{[N2]+ LN ]}, (9) 


then (8) becomes 


qi= Aw2S;..exp[— L5Aya{ {LNs +[N]} sec x dh][Ne2]. (10) 


(ii) The concentration of nitrogen atom produced by pre-dissociation process is 


given by 
go= Byigy Snexpl — [Biv sec x dh|LNel, (11) 
h 


where S:.. represents the number of mean total quanta available, By. the absorption 
coefficient which is independent of the wave length of the radiation, and gy the rate 
constant of the transition in (2) (b). 
§ 4. Formulae for the Equilibrium State 

Case I. Region in which the processes (1) and (4) are effective. 

The variations of [N2*] and [N] with the time are given by 


aN") _ 9,—al Neel, (12) 
ION) = alot Iel—eLNT, (13) 


where [e] represents the density of the electron per c.c., and «,, @ the rate constants 
of two body collision and dissociative recombination respectively. In an equilibrium 
state we have 


qi =m(NF, 


i.e., : 
Ay2S,.exp[—L5 Ay,{{INs1+IN]} sec x dh][No]J= (NF. (14) 


Case JI. Region in which processes (2) and (5) are effective. 
A variation of [N.] is given by 


AWE 
ANE) gy be PEM (15) 


where *, is the rate constant of three-body collision. In an equilibrium state it 
becomes 
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1.€., 


BraktySam expl—| By.[Nelsec x dh\Ns] = «LN PLM]. (16) 


§ 5. Calculations of Dissociation of Nitrogen 
Case I. 
we start to solve the equation (14). 

(14) becomes 


eol-1sA af LNA+ND sec x dh}=S0 NT (17) 
h Ay 29100 
We assume that 
dp=—{2[Ne]+LN ]} my gdh, (18) 
P={LNe]4+LN]} AT, (19) 


where p represents partial pressure, my an atomic mass of the nitrogen, g and k 
gravitational constant and Boltzmann’s constant respectively, and T the temperature. 
Then 


d 3 
ga CNA RT |=—@LNHVD my g. (20) 


If we take T as constant (20) becomes as follows, 


d 2 My Z_ oy ee Std 
Gr NAFLD =— CLINI DE = 2M IND Hg 2) 
When we put 
[N2]+[V]= a, (22) 
we obtain from (21) 
as. 1 
th ® = 10+ Ne]. 
Hence, after integration we get 
ete | . 
Q= zy-| [@an+- {fava lan, (23) 


where Q.. which represents the value of Q in the infinity is neglected compared with 
Q. : 
From (23) we get s a 

— \ Qah= — Qi + {[alah. (24) 


(i) Region of Small Dissociation 
In this region we have following approximate expressions 


Q = [Ne], (25) 
{2 an=| fA) an+|{Nydh~| Es) dh. (26) 


The limit of the dissociation that the approximation (26) holds is examined in Appendix 
II (a). Thus , 
— | ah=F QHy~—7 INA Hy. 27) 
h 
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(17) becomes using (22), (24) and (27) as follows, 


exp {—0.75 Aw2 HyLNe]} =p for” secy=). (28) 
(ii) Region of Large Dissociation 
In this region following approximations hold. 


Q =[N], (29) 
{iva dh<QHy. (30) 


The limit of the dissociation 
that the approximation (30) 
holds is examined in Appendix 
II (b). From (24) we get 


— {ean = ~QHy~—LN Wy. (31) 
Then (17) becomes as follows, 


eLNP 


Ags’ for secy=1. 


(32) 

If [NJ/LN.] is estimated, the 
values of [N.] or [N] which 
satisfies (28) or (32) are obtained: 
But in this method it is unable 
to compute the values of [N,] 
and [N] in the dissociation from 
about 20 percent to 50 percent. 
In the calculations of (28) and 
(32), S;.. is taken to be 2 x 10°/cm? 
/s by interpolation of the curve 
given by Nicolet [12]; «, is taken 
to be 10-'%cm'/s, 10-%cm*/s and 

10-"cm*/s as mentioned in Ap- 
pendix I, and as the value of Hy 
2x10%m, 3x10°cm and 5x 10%m 
are used. In the case that (28) 
and (32) cannot be used, the 
following method is used. That 
is, when [N2] and [NV] are known 
for the upper and lower limits 


of the dissociation we can know 


Fig. 1 Curves of [N,] vs. fQdh, [N2]/Q vs. fQdh and 
= [N2]/Q vs. [Ny] for x1=10-!8cm3/s, 10-%em3/s 
Ke dh by (25), (27) and (29), (31). and 10-%cm*/s when Hy =2X10%cm. 


immediately the values of Q and 


Accordingly we are able to plot the curve [N.1/Q vs. [No], [N21/Q vs. 


vs. 


[N22], Q and \ Qdh for various degrees of dissociation. 
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Qdh, 


ar 


or [N2] 


Qdh. By the interpolation of these curves we can easily obtain ANS values of 


mentioned above for Hy=2x10°cm are shown. 


After knowing the value of Q it is easily able to get the partial pressure., That 


is, 


p=QRT. 


In Fig. 1 (a), (b) and (c), curves 


(33) 


The values of [N2]/Q, [NV], [Ne], Q, {aan and p are shown in Table I for Hy =2 x 10%m 


v3 i Qdh (/em® column) [Ns] (/e.) 


=10-8 


%4=10-9 | y=10-8 


1.0 
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0.6 


0.5 
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ie 
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7.82 x 1016 
5.46 x 10% 
3.73 x 1016 


11.55 x 100"| 


5.42 x 10-17/4.15 x 10" 


2, 55% 10-3311.55,.x 104 


1.73 x 10-17/7.30 x 10” 


1,34 x 10-1” 
tA B75 x 16 


2.62 x 1016 


1.95 x 10% 


1.33 x 10% 


8.55 x 10" 


9.07 x 10-46/1.55 x 10° 


7.23 X 10-'67.00 « 10° 


5.74 « 10-16/3.10 « 109 
4.30 x 10-16/1.05 x 10° 


4.3110" 


2.80 x 10-16/2.40 « 108 


.60 x10" 


Table 1 


Q (/c.c.) 


%4= 10-19 | x; =10-18| x, =10-19 


PNG Gicic.) 


pb (dyne/cm?) 


%y= 10-8 |x = 10-19 


%4= 10-18 


y= 10-19 


5.42 X 10'4.15 x 1045.42 x 10% 


2.55 x 101.72 x 1042.83 x 104 


1.61 10"/9.12 x 10'|2.01 x 10" 


0 0) 


199 Sl0e2 


2.48 X 10-2 


1.72 X 10"2.83 x 10" 


1.82 «10"/4.02 x 10 


1.01 x 10"%6.55 x 10"/1.43 x 104 


7.00  10'/4.66 x 101.16 x 10"! 


7.88 x 10-3 


4.17 x 10-3 


1.32 x 10-7 
9.13: < 10-* 


1.95 x 1014.30 x 10” 


3.06 x 10-3 


6.56 x 10-3 


1.86. 104.64 x 10° 


oll 


5.10 10-3 


4.30 x 103.10 x 1018.60 x 101° 


2.35 X 10""11.75 x 10"|5.90 x 10° 


1.35 x 101.03 x 10"14.50 x 10" 


1.55 x 10°4.30 x 10” 


1.05 x 10'9)3.54 x 10” 


1.42 x 10-8 


8.1210 


3.96 x 10-8 


Phil s< ile 


7-21 KAO? 13.10 xX 108 


4.74 10-4 


2.06 x 10-3 


5.50 x 10° |5.31 x 10? 


2.75 x 10” 


4.25 x 10? |2.20 x 101 


2.41 10-4 


126 >< 105° 


1.60 x 109 |2.40 x 109 {1.60 x 10" 


2.16 «10° |1.54 x 10” 


1.12 10-4 


Tao xa Oes 


In this method we are unable to know the degree of dissociation at a certain height 
because the distribution [N.] or Q are not given, except the expressions (18) and (19). 
This difficulty is excluded using the data of measurement of the rocket flight in 
which the pressure and density are directly obtained. When we use (9), the observed 
pressure is given by 


Povs=1LSP=15QRT. (34) 


Using (34) we get the distribution of the nitrogen atom and molecule with height. 
Fig. 2 (a), (b) and (c) show these curves for H=2x10°cm, 3x10°cm and 5x10%cm 
respectively. The full lines show the distribution of [N] and the dotted lines show 
the distribution of [Nz] (only for (c)). It is found that Q and maximum [JN] for 
smaller Hy are greater than those for greater Hy. The curves of the distribution 
of [NV ] are alike one another and have maximum value at about 140km-150km 
altitudes, and dissociations begin at about 130km and proceed to completion at about 
150km-160km, 160km-180km and 200km-220km corresponding to #,=10-*°cm*/s 
10-“°cm?/s and 10="*cm*/s respectively. To determine an appropriate value of Hy, the 


78 


AN oy | sit i , eek (a) 


H,=2!0'om 


Fig. 2 


cS os 


V8 A Aq 
sl ee 


—+—+ 


(b) 
Hy=3xI0 on 


ip 10" 10 
CONCENTRATION OF PARTICLES 


Curves of distributions of Q (hair lines), 
[N] (full lines) and [No] (dotted lines, only 
for (a)) with height for x,=10-'%cm3/s, 
10-"cm$/s and 10-%cm3/s when Hy =2X 108 
cm (a), 3X10%em (b) and 5X10%cm. 


oxygen and dy, 
density. 


following method is used. As is 
shown in Fig. 2, above 200km it 
appears that [V]>[N.]. Then in that 
region the following relations are 
obtained. 
Povs= {LN ]+LOD}AT=3[O]FT, 
doys=[N] my +LO]m,, 


(35) 
(36) 


where m, is the atomic mass of the 


is the observed 
Thus 
T=[2 my+mMolPors/3kdovs 
=1.86x 107-4, (37) 
Obs 


2 Pons 
3kT* 
For the region from 190km to 220km 


T and W are shown in Table II. 


[V]= (38) 


' Table I 


h (km) | T (°K) | [N] (/c.c.) 


220 463 1.92 10° 
210 436 2.96 x 10° 
200 406 4.66 x 109 
190 372 7.20 x 109 


Now it is necessary that the value of 
[N] deduced in the calculation is 
consistent approximately with that in 
Table II. Hence it appears to be 


better to take it as Hy=2x10'cm. In conclusion, the distribution of [N] is given in 
Fig. 2 (a) and the maximum of [N] is of the order of 10" ec. 
Case II. 
In the case of the equation (16) the expressions (18). and (19) are also used. 
From (23) we obtain 


[tel = HyQ— 


, 


“Qah, 
h 


(39) 


and approximately [M]=3Q/2. When [N2]>[N] 
[;iNesdte=Hty0—| “Lean. 
h h 


Hence 


[twat nyo~d Hy tN. 


Substituting (40) in (16), we get 


(40) 


\ 
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L5eLNPLN>] 
By; LyS y 


When we give [N1J/[N2], [N2] which satisfies (41) is obtained. Rate constants By,, 
&y are taken to be 10-°cm:2, 


exp} — Bysy HylNs}} = (41) 


10-*°/s derived by Deb in one case 
and By2gyS:.=10-" derived by 
Bates is used in another case. *; 
is taken to be 1.5x107*(c.c.)?/s 
computed by Bates. S;.. is taken 
to be 3x10"/cm?/s from the 
curve derived by Nicolet. The 


HEIGHT IN KM 


Fig. 3 Curves of distribution of N with height in the results obtained using these 
lower atmosphere for By29¢yS:.=3X10-"™ (a) values are shown in Table III 


and By29S:.=10-" (b). and in Fig. 3. It is found that 
the value of [N] is almost the same magnitude through the levels from 130km to 
90km and very small compared with the value of Q. 


Table III 


| IN] 
Rs oe SSeS ee eee Q 
|Bwv28 xSse0= 10-"B wag ySi0=3X10- 


130(km)| 1.3 x10"(/c.c.) 4.0X103(/c.c.)) 3.3X10"(/c.c.) 
120 1.2x 10” 3.6 X 10° 7.610" 
110 1.1x 10" 3.4108 2.5102 
100 1.0x 10 2.8 x 108 ~ 1.11013 
90 7.8 10° 2.5 108 6.01018 
80 2.5 10° 1.2108 4.0 104 | 


§ 6. Discussions 

The results obtained in our method show that the dissociation of the nitrogen 
begins at the level of 120km-130km and proceeds to completion at the level from 
180km to 220km corresponding to the various values of rate constants used here. 
This is well consistent with the assumption of the dissociation by the Rocket Panel 
and Callmann [13]. The height of the beginning of the dissociation varies slightly 
for the conspicuous variation of the rate constant, while the height of completion of 
the dissociation varies greatly. Hence many constants used here must be studied in 
detail. Further the assumptions (18), (19) and (9) in the calculation should be minutely 
examined, though it is unlikely that the results vary remarkably. Thus the nitrogen 
above the level of about 200km appears to be almost in the atomic state, but this 
result is not yet conclusive. Bates [14] recently asserts, based on’ the spectroscopic 
experiment of the auroral spectrum, that the contribution of the molecular particle 
to the ionization in the F2 region is ynimportant, Judging from our results that may 
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hold. As the dissociation of the nitrogen may remarkably influence on the formation 
theory of the F1 and F2 regions, it hopes that more accurate computations are made. 
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Appendix I 


(a) The Range of Application of Dissociative Recombination and Pre-dissoeiation 
The effectiveness of the two dissociation processes is known by comparing 41 
arid q. From (10) and (11), q, and q, are written again for cosz=1, 


PDT or ime ew Av.l"Q dh\LNeI, 
h 
q:=S,.By.gexpl|—By J"uwa dhi\(NeI- 


In these expressions Aw., By», Zy, Sjo, So are given as follows, 


Avy2=10-"cm? Bates and Massey, 
By2=10-"cm? Deb, 
fy Sits Deb, 


Si0o =2X10fom'?/s  $...23%10"/em*/s Nicolet, 
or By2gvS;.= 107” Bates. 
At first, we adopt values given by Deb for By2, gy, then 


= 10"exp[-15 x 10-"{"@ dh). (42) 
2 h 


9 


- =1 when 15x10-"("Qah~I6. 
2 h 
16 


[5x 19-10". 


Hence, \,@ dh = 
h 
In Table I this is the value corresponding to the level [N.]>[N]. 


Thus HylNal=2{"@ dh=2 10°, 
h 


2x 10'8 


LNl= 3x0 


=10". 


When we use gyBy2S,..=10-" given by Bates, 


Ces asl fies 
fe = 3X10! expl-15x 10 ( Qdhi, (43) 
i =1 when 15x 10-"| “Qdh=85. 
2 h 
Hence [N2] = 56x10" =>0 x 10", 


19° 
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Consequently the dissociative recombination is more effective above the level where 
the total perticle is at least 5.6x10"/c.c.. 
(b) Comparison between Two-body Collision and Three-dody Collision 
The ratio of the recombination by three-body collision to that by two-body 
collision is given by 
roe (44) 
where x is the concentration of the total of the different kinds of particles, 0 is the 
ratio of the concentration of recombinig particle to . «, and «x; are given by several 
authors as follows, 
*,=10-'*(cm*/s) Herzberg [15], for oxygen atom, 
10-%—10-¥ Deb, for nitrogen atom produced by dissociative recombination, 
10724 Deb, for nitrogen atom produced by pre-dissociation, 
20>"5 Nicolet [16], for oxygen atom, 
and «,=10-**(cm*/s) Moses and Ta-You Wu, Deb, 
10-"*—10*% Penndorf, 
15x10-" Bates, 
5 UN Mitra [17]. 
We take «2 to be 10-**cm*/s and «,; to be 10-%—10-*°cm?/s and 107?!cm*/s for two 


cases noted above respectively. Then 
7=(10-"4~10-"5) n(1—8), 


or r=10-'n(1—8). 


Thus the recombination by two body collision is predominent above the level where 
the concentration of the total particle is at least 10'/c.c.. 


Appendix II 


(a) The Range of the Dissociation of Nitrogen for Application of (26). 


If we assume that 


=) 
— — H kT 
[Nino he N2]= flor Nye m » Hy.= 2myf? (45) 
was (iehy) 
oo oo od N. 7 
- then (cin =\" *—Frnye = Eh ey, (46) 
“0 0 


where a is constant, [N2} is the value at 4) and the upper limit is taken as infinitive 


instead of h=h +a, for convenience of computation. 


Then i [Ne]dh~LNob Hv2, (47) 
d ho 
and ~ QHy =2 Hwi Neh. (48) 


tis 2 
Hence QHy / J, uN] dh= Har" (49) 
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When a=100km and Hy2=10km (49) becomes 
r aye ~ ii" =20 
(Viwaah/ | CN dh~10, QHy J" LV 1dh=20. 
If h, is the level of ten percent dissociation, 


then \" wwelah / JS LN dh =5, Qn, Hy J [, (Nah =10. 


(b) The Range of the Dissociation of Nitrogen for’ Application of (30). 
If we assume 


h-ho 
[NsJ=LNsbe |”? ¥ 0) 
then \ LNel=HvLNed (51) 
0 


If degree of dissociation is fifty percent at fr, 


HyQ=4Hyx{LNeh; (52) 
and | [[Neldh/Hy Q = (53) 


As it is estimated that [N.] decreases more rapidly than that expressed by (50), the 
ratio in (53) becomes smaller. Hence it seems that above the level of dissociation of 
fifty percent (30) approximately holds. 
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Propagation of the Cosmic Radiation through Interstellar Space 
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Department of Physics, Department of Physics, 
Osaka City Univervsity Osaka University 
Abstract 


Various theories on the origin of the cosmic radiation are 
critically surveyed. Expected phenomena resulting from the galactic 
origin theory in which charged cosmic rays are trapped by the glactic 
magnetic field are discussed with special regard to the secondaries 
due to the collisions with interstellar hydrogen. Electrons through 
the decay process z>p—¢e are then expected to amount 3~9% of the . 
total intensity, depending upon the magnetic field strength in our 
galaxy. An alternative hypothesis that the cosmic radiation propagates 
in a straight way is also examined. It is proposed that this may be 
tested by observing primary photons, resulting from the decay process 
m°->2r, which could be several percent of the total radiation at low 
latitudes. 

1. Introduction 

It has been thought plausible that the cosmic radiation is stored in a finite 
region of the universe, since otherwise the amount of energy carried by the cosmic 
radiation is far larger than that by the electromagnetic radiation that might be a 
source of the former energy. This has, on one hand, led Richmyer and Teller (1) to 
their speculation that the cosmic radiation is confined only in the solar system, and, 
on the other hand, Fermi (2) and Unsdéld (3) to their theories of its galactic origin,*) 
where the cosmic radiation is kept inside due to the wandering magnetic field sug- 
gested by Alfvén (4). 

Although these theories are ingeniously elaborated, a number of difficulties have 
been poined out. The solar origin does hardly harmonize with the unappreciable 
diurnal variation of very energetic cosmic rays which can not so much curl in the 
magnetic field as to be made isotropic (5). The galactic origin, though interesting in 
view of their acceleration mechanism, meets an objection, when experiments on heavy 
primaries are taken into consideration. A negligible amount of Li, Be and B in the 
primary cosmic radiation observed by Bradt and Peters (6) led them to conclude the 
small thickness (S2gcm~-*) of matter traversed by the cosmic radiation from its 
origin to the top of the atmosphere, while the galactic origin theory requires the 


*) By the galactic origin we mean that the cosmic radiation is created in our galaxy and stored 
in it, if not specially noticed, 
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thickness of traversed matter as large as the absorption mean free path of cosmic 
rays (~40gcm~’). 

The scarcity of Li, Be and B is, however, opposed by a later experiment by 
Daniton et al (7) who observed as many Li, Be and B as C, N and O. This allows 
the thickness of traversed matter to be larger, but its upper limit is set by the 
existence of heavier nuclei whose absorption mean free path is about 10gcm~*. Thus 
the contradiction with the galactic origin hypothesis is yet unsolved, irrespective to 
whether the abundance of Li, Be and B are appreciable or not. 

Although these experimental results are unfavourable to the galactic origin, it 
seems worth while to seek for further evidences for or against this hypothesis. This 
is our principal aim of this paper, a preliminary account being published in a previous 
short note (8). 

2. Diffusion of the Cosmic Radiation in the Galaxy 

According to Fermi (2) and Unsdéld (3), cosmic ray particles are scattered by 
magnetic clouds whose linear dimension, LZ, is of the order of 10'%cm (Fermi) or 10*°cm 
Unsoéld. If the cosmic radiation is created at r=0 with unit intensity, its density at 
y after time ¢ is given by (9) 

o(7, t) =(274?)-*/* exp{—(7?/2 2°)—t/T}, (2.1) 
where the boundary condition is natural and different from that assumed by Cocconi 
(5). A is the diffusion length and is related with the number of collisions with 
magnetic clouds, N, as 

AM=2 NEA 3 ]2ctL/s. (2.2) 
T is the mean life time due to the colliisons with interstellar matter, mainly with 
hydrogen of density 1 atom cm™=*. 

With the collision cross section for protons ~5x10-2%cm?, 


T=7x10"sec =2 x 10’ y.» (2.3) 


For an instanteneous source, such as suggested by Sekido ef al (10) the density 
reaches the maximum at 


t= TWV147/3eLT—1}. (2.4) 
In Sekid’s example, the coimic radiation would have been originated 1.2 10"'sec ago 
at 3x10"'cm apart from the earth (at the Grab nebula). Then the cosmic radiation 
would be far before its mean age and its intensity would be increasing with rate 
3c/4L, corresponding to 0.2% y-! for L=10%cm. Thus the contradiction between 
Sekido’s (10) and Cornell’s (5) (11) results should not be due to the difference in the 
periods of observation if Sekido’s conjecture were approved. 


The density distribution (2.1) results in anisotropy of the cosmic ray intensity 


—(L/3) p/p =r/2ct. (2.5) 


*) The absorption mean free path may be only slightly different from the collision mean free 


path in such diffuse matter, where pions takin i i i 
; g main part in N-cascade i 
undergo decays before collisions, en 
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This leads in Sekido’s example to the anisotropy of about 3/8, which is, of course, 
covered by a large background. If the observed anisotropy of about 10% were true, 
the Crab nebula should contribute to the total intensity by nearly 22%. Since the 
contribution is far before its maximum at the present, its contribution would be 
unreasonably large in future. If the explosions of supernovas such as the Crab nebula 
were to account for the whole intensity of cosmic rays (12), the explosions of every 
300 y would result in too large intensity to be consistent with our experiences, 
provided Sekido’s work were to be acceptable. 

If the cosmic radiation is continuously produced with a constant intensity, that 
produced at T before gives the largest contribution. For such a source locating at 
the center of our galaxy, the anisotropy is expected to be about 0.07% which may be 
too small to be detectable by present techniques. It may be noteworthy that the 
maximum of the cosmic ray intensity expected is just in an opposite side of the 
observed one (13), though the observation on siderial time variations is not conclusive 
yet. - 

There could be another cause of anisotropy, if the source is at the center. 
Since the solar system flies away from the center with the velocity of about 
3 x10°em sec™', the intensity at the opposite side could be larger than that at the side 
of the center by about 0.02%. This is again too small to be observed. As seen above, 
the cosmic radiation should distribute practically isotropic in the galaxy, independent 
of its positions of sources, so that the uniform density may be assumed in the follow- 
ing qualitative discussions. Then the density distribution (2.1) contains only a time 
dependent factor exp (—#/T). 

Now we consider the secondaries produced by the collisions of primaries with 
interstellar hydrogen. The density of the secondaries is given by 


Ppp en es), (2.6) 


where T: represents the mean life of secondaries and f the fractional intensity of 
secondaries produced per collision with the same energy as the incident primary. For 
large ¢ when the equilibrium between primaries and secondaries is reached, the ratio 
of secondaries to primaries tends to 

Jfe/ Tf. (2.6') 


This is employed by Bradt and Peters (1) in estimating the contribution of heavy 
primaries. 

So far we have assumed that the cosmic radiation originates solely in the 
galaxy. The case may be plausible where a part of the cosmic radiation comes from 
out of our galaxy and is trapped due to the galactic magnetic field. Assuming that 
the intensity outside the galaxy is J and a fraction, g, of the cosmic radiation striking 
the galactic surface can penetrate, the density inside the galaxy follows the equation 


d__(1,&8 S 
oo -(4+ Spel )ot2n ely (2.7) 
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where S and V are the surface area and volume of the galaxy, respectively, and a 
term cSg’/2V expresses leakage through the galactic surface. (2.7) is solved as 


9=2n(S/V) ghTi(l—e-/"1) (2.8) 
with 
1/T,=1/T+¢eSg'/2V. (2.9) 
For large ¢ this tends to 
p=2n(S/V) ghT:=s Zoo, (2.10) 


where p)=471)/c is the density outside and s=cT,S/2V may be called the storage 
factor after Unsdéld (3). On account of (2.9) the storage factor is given by 


$=8)/(1+S) 8"), (2.11) 


where s,=cTS/2V is the storage factor without leacage. Since s, is of the order of 
10°, the density inside is much larger than that outside, provided sog’ <1. Then sources 
outside contribute to the observed intensity by an unappreciable amount and the 
supernova origin theory (12) would have to be reexamined for being reconciled with 
the galactic magnetic field hypothesis. 

3. Secondaries Produced by the Collisions with Interstellar Matter 

The collisions of heavy nuclei with interstellar matter have been discussed by 
Bradt and Peters (1). Fan (9) argued that this process could account for the observed 
abandunce of heavy nuclei with the acceleration mechanism proposed by Fermi (2). 
The energy spectra they obtained should be steeper, the heavier the nuclei are, in 
contradiction with recent observations (14) (15). This is another difficulty against the 
galactic origin. Further difficlties are found if other possible secondaries are taken 
into consideration, as will be shown in what follows. 

Since most of the primary cosmic radiation are protons, it may be sufficient for 
our purpose to deal with the collisions of protons with interstellar hydrogen. As the 
order of magnitude theory, the intensities of secondaries may be estimated in con- 
sidering the minimum energy of primaries responsible to the production of the seconda- _ 
ries under consideration, as has been done in our previous note (8). In this method of 
estimate was assumed the energy spectrum of the secondaries not appreciably different 
from that primaries. This is really a good approximation at high energies and for 
the processes we consider. For the purpose of definiteness, however, we shall calculate 
the spectra as well as intensities of secondaries in reference to Fermi’s theory of 
multiple meson production (16), though the theory is only tentative and other theories 
of meson production lead us to similar conclusions. 

For later purposes we are mainly concerned with thin interstellar matter and 
the cascade development of protons is discarded.) Then the energy spectrum of 
protons changes so little that we may assume the spectrum everwhere to be the same 
as that observed at the top of the atmosphere. This is given at high energies in a 
good approximation by (15) 


p(E)dE=p,E-"dE, (3.2) 


*) See also footnote at page 84 
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where £ is the total energy of a proton and @=2.0~2.5 gradually increasing with 
energy. The use of such a power spectrum is suggested by analytical convenience. 
Secondaries produced by proton-proton collisions are protons, neutrons and their 
anti-particles, if exist, and also various mesons. Among them only protons, anti-pro- 
tons, electrons, positrons and photons are stable, the latter three being dacay products. 

(a) Anti-protons. The search for anti-protons in the primary cosmic radiation 
would provide an interesting test for or against their existence, if their intensity were 
high enough. In the galactic origin theory, their intensity is expected to be f times of 
the intensity of protons, because T,= Tin (2.6). There are neither reliable experiment 
nor trustable theory for the magnitude of f, however. We can only set an upper limit 
of the intensity based on a primitive consideration. 

The threshold energy for producting an anti-proton is 8 Gev. Near the threshold, 
however, the volume of momentum space of the final state is so small that the cross 
section is expected to be very small. The cross section may reach as large as the 
geometrical one only above 100 Gev. Thus their intensity is at most of the order of 
1% of the proton intensity. Taking into account the competition with the production 
of mesons, the actual intensity may be a fraction of percent. Hence the existence of 
anti-protons could hardly be examined. 

(b) Electrons (positrons inclusive). Neutrons are possible to contribute to 
electrons through their §-decays, but the energy imparted to an electron is so small 
that they are neglected in comparison with other processes. The largest contribution 
comes from charged pions through muons, while other unstable particles may 
contribute to small percentage. 

For comparsion with the latest experiment at geomagnetic latitude 55° (17), we 
estimate the minimum energy of primary protons contributing to electrons whose cut 
off energy is 1.7 Gev. Since an electron takes on the average 1/3 of the energy of 
a parent moun, the minimum energy of muons under consideration may be taken as 
about 5 Gev. Those muons are the decay products of pions with energy greater than 
6 Gev, because the average energy of a muon is about 0.8 times an energy of a parent 
pion. If 2 pions with average energy E, are produced by a proton of energy E, 
there holds a relationship 

27(r-l=nE,, (3.2) 
provided pions are produced istropically in the center of mass system. 7 is the Lorentz 
factor that transforms the laboratory system into the center of mass system and is 
expressed as 

y=V (E+1)/2, (3.3) 
where the total energy E is measured in unit of proton rest energy. Then we obtain 
E>10 Gev for n=1, E>17 Gev for n=2 and so forth. For E=10~20 Gev the average 
number of pions produced may be two. Hence we can roughly set the minimum 


energy of protons for producing observable electrons as 


E=20Gev, (x—p-@), (3.4) 
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taking into account recoils of colliding protons. This gives to fraction f in (2.6) 
Ff =0.05, (x—p—@) (3.5) 


for B= 2.0 in (3.1). Thus the intensity of electrons would be appreciable, if protons 
were to traverse thick enough intersteller matter. 

A more quantitative estimate is carried out in reference to the Fermi’s theory 
of meson production and to a power spectrum like (3.1). The differential spectrum 


of electrons is obtained as 
hy (E., t) dE, =0.080 py (1— e*/) E78 9/3 dE... (3.6) 


E, is the energy of an electron in unit of proton rest energy. The intensity of the 
electrons is obtained by integrating (3.6) above a magnetic cut-off energy E,° as 


H,(E.°, t)=(0.24/(4 B—5)) po (1—e-/7) EGP 8, (3.7) 
Thus the ratio to the intensity of protons is given for E£,°=1.8 as 


(i (8=2.0), 
(0.017 (@=2.5). 


(3.8) 


These figures are in rough agreement with (35). 

If the electrons thus produced are trapped by our galaxy with mean life 
time 7; their intensity is given by (2.6) with (3.8). The mean life time T, is subjected 
to the interactions of electrons with galactic electromagnetic wave, as discussed by 
Feenberg and Primakoff (8) and elaborated by Donahue (19). According to their 
investigations the most important process for the energy loss of electrons is the 
Compton collisions with galactic photons. 

Since both the cross section and energy loss for the Compton process decrease 
above the electron rest energy mc’, only the nonrelativistic limit in the rest system 
of an electron may be considered. Then the cross section is that for the Thomson 
scattering and the energy loss is °e, where € and mc? is the energies of colliding 
photon and electron in the observed system. With the energy spectrum of photons 
n(€)de, the energy loss of an electron is approximately evaluated as 


dE,/dt=—co,{ SOC i (e) &e de, 
“0 


o,= (87/3)(e?/me")*. j (3.9) 


On account of that n(e)ede is the differential energy density spectrum, the energy loss 


(3.9) is proportional to the energy density of photons with energies below mc?/é 
Hence 


adE.,/dt = —co,€* W(me2/&) = —2 x 10-42 W(me?/é) ev sec", 
W (me?/é)= es : 


‘ n(e)e de. (3.10) 
W(mce?/é) may be equated to the energy density of galactic photons. This is 


about 0.3 ev/cc for visible light and the order of 10-‘ev/cc for radio frequency 
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radiation. Hence only the visible light is taken into consideration. 


The life time for an electron of energy € mc? losing its energy down to &mc? 
is then obtained as 


T.=5 x 10"(me?/ W) (1—€/&)/€ sec. (3.11) 
For &)mc?=10 Gev and £me?=5 Gev 
T.=4x 10" sec. (3.11/) 


This is larger than the mean life time of protons. 

We may, therefore, expect a considerable intensity of electrons with energies 
below 10 Gey in the primary cosmic radiation. The differential energy spectum of 
such electrons that are produced by the collisions of protons with interstellar hydrogen 
and suffer the collisions with galactic photons can be obtained by solving a diffusion 


equation 
: OE, a Abe h(E,, t) 


= k,— *, 


as t) 4+S(E, 0). (3.12) 


k.E.? is the energy loss given by (3.10) ~_ 
k.=co.W/(me*)*=2 x 10-*ev—'sec™'. (3.10’) 


The continuous energy loss adopted above may be permitted on account of that several 
Compton collisions take place in the mean life time of-an electron and the energy loss 
per collision is much smaller than E, in our case. The source of electrons, S(E., ¢), 
is taken from (3.6) in an equilibrium state as 
S(E., 1)=(So/T)E", 
u=2(28—1)/3. é (3.6’) 


The diffusion equation is soved as 
WE. t)=S,E.-*-1/k.(u—1). (3.13) 


The intensity of electrons is obtained by integrating (3.13) for E, larger than the 
cut-off energy, E£,°. 

HE.) =S.Eo-"/k,u(u—l). (314) 
The ratio to the intensity of protons is about 0.09 for #,°= 1.8 and #= 2.0, which is 


large enough to be detectable. 

If there were strong magnetic fields in our galaxy, the radiation loss of electrons 
due to curling in the magnetic field would not be negligible either.” The radiation 
loss of an electron with energy E, in a magnetic field of strength H (in gauss) is 


given by 


oie =—1 x 10-°H2é" evsec"!. (3.15) 


Since this process has the same energy dependence as the energy loss by the Compton 
effect, one has only to add a term ky to k,. Comparing (3.15) with (3.10), and (3.10), 


we have 


*) The authors are grateful to Professor Rossi for callidg our attention to this process. 


‘ 
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kg=1x10°H?/(me*)?=4 x 10" H?ev-" sec“. (3.15) 


Hence for H~7 10-* gauss the radiation loss is comparable with the loss by the 
Compton effect. For the maximum possible magnetic field strength, 10-° gauss, ky= 
4x10-2° results in the reduction of the intensity of electrons by a factor of three in 
comparison with the above estimate. In this case electrons should be barely dectable 
in the primary cosmic radiation. 

The absence of electrons in the primary cosmic radiation is, therefore, found 
to be a clue to testify the galactic origin. There might, however, be ambiguities in 
the above estimate, although all but the density of galactic photons do not seem to 
depress the ratio of the electron intensity to the proton intensity below 1%. We have 
assumed the photon density as uniform throughout the galaxy, but it is merely a 
rough approximation. If the electrons spend more: time in strong magnetic fields 
where the density of photons might be larger than in other parts, the intensity of 
eletrons would be smaller than that given above and our argument would have to be 
reviced. In this case, however, the evidence on heavy primaries remains to be against 
the galactic origin. 

(c) Photons. In the experiment referred (17) photons of energies greater than 
0.73 Gev are tried to observe. The energy of a neutral pion which produce photons 
of this energy is about 1.5 Gev. This pion can be produced by a proton of energy 
above 3 Gev. Taking into account the recoil of nucleons, the minimum energy of the 
protons responsible to the observable photons may be assumed as 4 Gev. Hence their 
intensity at 55° is about 1/2 of the total proton intensity. Around this energy most of 
pions are produced singly and 1/3 of them may be neutral. Thus we obtain the 
fraction of observable photons per proton-proton collision as 


F-=(2/3) x (1/2)~0.3. (3.16) 


More quantitative estimate is made in reference to Fermi’s theory as in (b). 
The differential spectrum of photons is evaluated as 


g(E;) dE, =0.35 p)(1—e-/?) E22 8-3 (3.17) 


Integrating (3.16) above the cut-off energy E,° that depends upon experimental 
conditions, we obtain the intensity of photons 


L 
GE )=z fo po(L—e-!/?) Eye- 48-3973, (3.18) 


For E’=0.8 


( (A=2.0), 
(3.19) 


0.3  (@=2.5). 


In order to account for the intensity of the photons which are in equilibrium 
with protons, we must modify (2.6), since photons are not trapped by the magnetic 


field. The mean life time for photons is now estimated on account of escaping of 
photons out of our galaxy as 
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T2, ,=2V/cS~10" sec. (3.20) 
Hence the fractional intensity of photons is given by 


FAT3,,/T)=5x 104, (3.21) 
which is too small to be detectable. 

From the above discussions we could expect the detectable intensity of electrons 
against observations if there were such a galactic magnetic field trapping charged 
particles that were not too strong. Together with other evidences mentioned in 
Section 1, doubt may be thrown upon the existence of such a magnetic field as well 
as the theories of Fermi and Unsdéld. It seems now necessary to examine the 
hypothesis that the cosmic radiation propagates in a straight way through the galaxy. 
How one can test the latter presumption will be considered below. 


4. Possible Observable Effects for the Cosmic Radiation Going in a Straight Way 
through the Galaxy 

Under this assumption we may expect the anisotropy of the cosmic radiation, 
since the earth locates near the edge of our galaxy and the amount of matter 
traversed is different depending upon the incident direction. At siderial time 18h the 
cosmic radiation traverses the thickest matter of about 0.5g/cm*. This results in the 
absorption of several percent for heavy neclei. In practical cases the finite solid angle 
of an apparatus makes this percentage as small as the order of 1%. 

As for secondaries such as electrons and photons, (2.6) is applied only at small 
-t, so that their intensity is given by 


FH/A (4.1) 


where x is the thickness traversed in gcm~* and A is the collision mean free path 
for primaries. For protons traversing through the center of the galaxy, x/A=0.5/30= 
0.02. Now we can estimate the intensity of various secondaries as in Section 3. 

(a) Anti--protons. f is only the order of 1% so that this is out of question. 

(b) Electrons. This case is nothing to do with the absorption as discussed in 
Section 3. Thus electrons are expected to be practically unobservable either. 

(c). Photons. Now the large f, in (3.19) is of considerable importance. At 55° 
where the cut-off energy for protons is low, the intensity of photons is at most a 
fraction of percent, in cosistent with the observation (17). At low latitudes, hewever, 
the intensity of charged primaries decreases, whereas the intensity of photons remains 
unchanged. Then the relative intensity of photons increases with decreasing latitude 
and could be several percent of the total primary intensity at the equator. This is 
certainly observable by the present evperimental techniques. The directional 
dependence or the diurnal variation of photon intensity may also be observable. 

5. Concluding Remarks 

We have considered the various effects of collisions of the cosmic radiation with 
interstellar matter. On the basis of the galactic origin theory that charged particles 
are stored in our galaxy, a difficulty arises on account for the negligible intensity of 
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electrons which should be produced by the collisions with interstellar matter and 
trapped by the galactic field, though this evidence has to be taken with reservation. 
This and other difficulties regarding with heavy primaries may be unfavourable to 
the presence of the galactic magnetic field that traps charged cosmic rays for long 
periods, though an astrophysical evidence is reported for supporting its presence (20). 

An alternative hypothesis that the cosmic radiation is not stored in our galaxy 
but comes mainly from other galaxies seems rather plausible. In order to see whether 
this hypothesis is acceptable or not, the observation of photons in primary cosmic 
rays at low lattitudes is suggested. 

If this is approved, we shall have to turn towards the nonequilibrium theory of 
the cosmic radiation, as often presumed by Lemaitre (21) and others, because it takes 
as long time as the age of the universe that the cosmic radiation travels from the 
farthest galaxies to the earth. Then the composition and the energy spectrum of the 
cosmic radiation might be correlated with the origin of the universe. The difference 
in the compositions of cosmic rays and existing elements may be attributed to the 
acceleration of the former due to the gravitation effect that should be of greatest 
importance at the origin of the universe and is more effective for heavier elements. 
The electromagnetic acceleration is proportional to Z/A in ordinary conditions and 
hardly account for more abundant high Z elements in the cosmic radiation than in 
whole elements and for the unexpectedly large energy content of the cosmic radiation. 
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The Vertical Distribution of Electrical Conductivity in 
the Upper Atmosphere* 


By Hiroshi MAEDA 


(Geophysical Institute, Kyoto University) 
Abstract 


Using available materials of recent rocket measurements, a 
calculation is made to obtain the atmospheric and ionospheric models 
of the upper atmosphere between 60km and 400km. On the basis of 
these models the vertical distribution of electrical conductivity and its 
height-integrated values for three ionospheric regions F, E and D are 
obtained. The results show that the conductivity of the E-region, 
especially near 100km, is most predominant, so that this region is 
most probably effective to the electrical currents producing daily 
magnetic variations. 

1. Introduction 

Recently, from four sources [1] anisotropic conductivity of the upper atmosphere 
has been proposed simultaneously. This new concept of conductivity plays probably 
an important role in the electromagnetic interpretation of the problems relating to 
the upper atmosphere. Therefore, it needs to try a detailed calculation to obtain 
the vertical distribution of conductivity in the upper atmosphere. For this purpose, 
however, the atmospheric and ionospheric models must be determined. 

Since 1946, a relatively great number of sounding rockets have been fired in order 
to obtain direct information on pressure, density and temperature of the upper 
atmosphere. From these results the Rocket Panel [2] has reported a combined result, 
and given two tables of physical properties of the atmosphere up to 220km. Recently 
H.K. Kallmann [3] has studied theoretically the physical state of the atmosphere 
between 80km and 250km, and given a table which is in fair agreement with the 
results of rocket measurements. 

Very recently, the results of direct measurements of the electron density in the 
ionosphere by rocket-borne instruments have been reported by W.W. Berning [4], and 
J.R. Lien, R.J. Marcou, J.C. Ulwick, J. Aarons, and D.R. McMorrow [5]. 

In the present paper, the author shows the results of calculation of the 
vertical distribution of electrical conductivity and its height-integrated values, using 


atmospheric and ionospheric models deduced from the rocket measurements mentioned 
above. 


pe a a eee Oe es ee 
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2. Atmospheric model 
(i) Temperature (T) 
The vertical distribution of temperature up to 80km follows the result 


adopted by the Rocket Panel [2], and above 80km the following linear distribution 
assumed by Kallmann [3] is used: 


T=T)+r(h—In) (1) 
where 
hy =datum level=80km 
T,=temperature at the datum level=190°K 
ry. =lapse rate=2.5°K per km. 
The distribution is shown in Fig. 1. 
(ii) Molecular weight (M) | 
_ According to Kallmann [3], if the following assumptions are made: 
, a) Oxygen dissociates in the region between 80 and ~120km, 
b) Nitrogen dissociates in the region between ~130 and 250km, 
c) The temperature and pressure at the datum level (80km) are known and the 
lapse rate is 2.5°K per km, 
d) The fractional change in concentration is of the form 


dC(O, or Nz] _ 
C[O, or Nz] me Kidh 


where C is the concentration of molecular oxygen or nitrogen at any height, 
and K; is a concentration coefficient equal approximately to 1.93x10-%cm-' 
for O, and 0.98x10-°cm~! for N2, 
then the molecular weight as a function of temperature and of height is obtained 
from the equations 


dM, _[ M,—238 1 Mg -K7| : 
an oe 1 r+ K; (2) 
dM, M,—14.4 Mog | 

o | r | R +7r—K, (3) 


where M, and M, are the mean molecular weight of air in the regions where oxygen 
and nitrogen dissociates, respectively, and g is the gravity. Also in our calculation, 
these equations are used, and the result of calculation is shown in Fig. 1. 

(iii) Gravity (g) 

The Distribution of gravity with height is calculated from the inverse square 


law 
ee eek es 4) 
: e=0i( ath ) ( 
where a=radius of the earth=6372 km 
Z= gravity at the earth’s surface =980cm.sec.~* 


(iv) Scale height (H) 
Using above quantities, the scale height can be calculated from the relation 
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weight, and scale height. 


where & is Boltzmann’s constant. 


(viii) Geomagnetic total intensity (F) 
The geomagnetic total intensity is calculated from the inverse cube law 


F=F( 50) 


ath (9) 


where Fy)= geomagnetic total intensity 
at the earth’s surface=0.5 gauss. 
3. TIonospheric model 

(i) Electron density (n,) 

The results of direct measure- 
ments of the electron density in the 
ionosphere mentioned above are fairly 
complicated, so that, based on these 
results, we have assumed a model 
distribution of the electron density at 
daytime shown in F ig. -2. 

(ii) Jon density (7;) 

Recently preliminary rather 
succesful attempts to determine the 
ionization in the E-region by means of 
a probe technique are described by 


H. Maepa 


where # is gas constant. The distribu- 
tion of H is shown in Fig. 1. 

(v) Pressure (P) 

The pressure is obtained by 
numerical integration of the equation 


(6) 


(vi) Density (e) 
The density is obtained from the 
equation of state 


_PM 


P= RT (7) 


(vii) Number density (n) 
The number density shown in 
Fig. 2 is calculated from the relation 


(8) 


10 10 10) to 
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10° 
NEGATIVE ION TO ELECTRON RATIO 


Fig. 2—Altitude vs gas density, electron and ion 


(negative and positive) densities, and 
negative ion to electron ratio. 
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Gunner Hok, N.W. Spencer and W.G. Dow [6], and they found that a positive-ion 
density is about ten times larger than the electron density. This result is very 
interesting. But, for the purpose of our calculation the height-range of measurements 
is too small, so that we use a theoretical estimation of 4 (the ratio of the concentra- 
tion of negative ions to that of electrons) by D.R. Bates and H.S.W. Massey [7], until 
above direct measurements by rocket are extended over all altitudes in the ionos- 
phere. 

From theoretical considerations, Bates @nd Massey have estimated the value of 
4 in the different ionospheric layers (F2-, Fl-, E-, and D-layers). On the basis of 
these values, we have assumed the distribution of 4 shown in Fig. 2. Using this, the 
positive and negative ion densities are calculated under the assumption that atmos- 
pheric air is electrically neutral as a whole, and these distributions are also shown 
in Fig. 2. 
4. Collision frequency 

S. Chapman and T.G. Cowling [8] and T.G. Cowling [9] have given formulae to 
determine collision frequencies, and recently M. Nicolet [10], using the formulae 
derived by above authors, has calculated and discussed the electron collision frequency 
of all ionospheric regions. We follow also these workers. 

(a) Collision between neutral and charged particles 

In the velocity distribution method, the effective collision interval ti. of mole- 
cular mass 72; and mass m; is 


mym2(1+ 22) 
(7,2, + 222) kT (10) 


Fis — 


where D is the ordinary coefficient of mutual diffusion of the two kinds of molecules 
in the absence of other gases, 2, and m, are the particle number densities, k is Boltz- 
mann’s constant, and T is the absolute temperature. To a certain degree of appro- 
ximation, D is, for rigid elastic sphere, 


16 (my + Nz) 012" TIN, M2 


D 


where oz is the collision distance. 
From (10) and (11), the collision frequency v;, for a molecule of type 1 with 
molecules of type 2 is 
_16 mm+nm, of (+222) kT |r 
catia SS OT as [ 20m Mz ey 


(i) Collision frequency of an electron with neutral particles (Yen) 
For collisions of an electron of mass m, with neutral particles of mass m with 


number density x, (12) becomes 


_4 o( BRT \i 
Ven = ZNT en (= ) ' (13) 


because 
M=M> MFM, 
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Thus, the collision frequency of an electron with neutral particles is determined if we 
can find the effective collision distance. Several workers [9, 10, 11, 12] have baa 
v., using the value of o,, estimated theoretically or experimentally, but sip not easy 
to obtain an exact value for the collision distance of slow electrons in the ionosphere. 
We have, therefore, adopted the following relation to obtain v,, as a mean state: 


Ven = 5X 107 nV T (14) 
(ii) Collision frequency of an zonqwith neutral particles (v;,) 
Assuming that a mass of ions, negative and positive, is equal to a mass of 
neutral particle, i.e. m.=m,=m, (12) becomes 


“oe 1 
Vin= pmo Ar) (15) 


Following Cowling [9], we have used the value 


_e/ 300 16 
O1n= 59X10 ( er) (16) 


which is derived from measurements of the mobility of N.* ion in molecular nitrogen 
by A.M. Tyndall [13], and obtained 


Vin = 3.35 X 10-§ (17) 


(b) Collision between charged particles 
According to Chapman and Cowling, the effective collision interval for charged 


particles of mass m, and mz is still to be found by (10), but now diffusion coefficient 
D for charged particles is 


ee Oe eee ‘(ry es 18 
D=%5 rea Mz | e A,(2) (18) 


where e denotes the electric charge as used in the following equation for the force 
P between pairs of particles having charge e, and e: 


pa Lie 
ri 


and A,(2) denotes a slowly varying function of T and N,. 


(i) Collision frequency of an electron with ions (Yes) 
If we use the condition 


=m ;=m> V1 P11 2) 
then the collision frequency of an electron with ions of number density 2; is given by 


4 me? 
3 (22m,k*T*)3 


Veg = 


A,(2)n; (19) 


where the correction term A,(2) is given by Chapman and Cowling as follows: 


Aj (2) =log.(1+v9;°) 
where 
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with d=mean distance between pair of neighbouring molecules. 
If we define d as being equal to Debye’s distance, following Nicolet [10], 


my ee 
om ( 0) 


and if a | oe tht 
e*(2n,)8 
then 
a 4. (kT \i ; 
A,@)=2 log.= “weit ; ) (20) 


Using this correction term, the collision frequency of an electron with ions is given by 


3 


phys 
vei = E + 8.36 lossy fs Apa (21) 
z 


(ii) Collision frequency of a negatwe ion with positive ions (v_,), and of a 
positive ion with negative ions (v,-) 
Assuming that relation (20) holds in these cases and that m_=m,=m for simpli- 


city, then the collision frequencies v_, and y,_ are given by 


3.06 « 10-4 


be oe (22) 
-14 
1 = ven (23) 


Using these, the effective colli- 
sion frequencies y,, v. and v, are 
obtained by 


Ve=VentVeuotVex (24) 
as ee ae (25) 
Vy =VantVy- (26) 
and the distributions are shown in 
Fig. 3. 
5. Electrical conductivity Ee oilition Hpabuincy Why 8S alia 
Using the quantities prepared Fig. 3—Altitude vs collision frequencies. 


above, we can calculate the conductivity as below. 

(i) According to M. Hirono [14] and D.F. Martyn’s [15] expressions, if H may 
be resolved into components EH, and E, respectively parallel and perpendicular to H 
and hk is a unit vector parallel to H, then the current density is expressed as 


Qe obo tok +o(hx EF) (27) 

where 
a ee nN nN, 28 
cake Fas mv." mae) cee 
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_ eF er m3 
OO an 7 1. “mm 


10 16 10 10" 10 
CONDUCTIVITY (e.mu.) 


Fig. 4—Altitude vs electrical conductivities. 


At the magnetic equator, if an east- 
ward electric field is impressed, then 
the Hall current due to this is vertical, 
and so makes charges accumulate 
above and below the conducting layer. 
If they build up until the resulting 
vertical polarization field stops further 
flow of the Hall current, the resultant 
eastward conductivity is increased to 
ox( =a, +2), Fig. 4 shows the varia- 
tion of 0, 01, 2, and o3; with height. 

(ii) Let x-, y-, z- axes be taken 
to magnetic south, east and upward, 
respectively, then the current density 
J is resolved as follows [15]: 


J.=(acos’¢ + o,sin’$)E,,+ o.Ey sin ¢+(o—o1)E.cos ¢ sing (31) 
Jy=—E,sing+o,Ey+o:E,cos¢ (32) 
J.=(00—01)E, cos ¢ sin $—o,E,cos ¢+(o)sin?¢+o; cos*¢)E, (33) 
where ¢ is the magnetic dip angle. 
Assuming that J,=0 at any height in the ionosphere, then 
_ 62,Eycos $—(a)—01)E,cos ¢ sing 
E, a osin*¢ + a\Ccos*¢ (34) 
so that, we have 
ie = Tuglin = OxpyEy (35) 
Jy=—OyE, +OyyEy (36) 
where 
ee Foo 
*% dysin®¢-+o,cos*¢ (37) 
- doo2Singd 
709 G SiO 0:C05H = 
__ 9998in?¢ + (¢,2+ 02") cos’ 
ov = a sin*¢ +0,cos?¢ (39) 
(iii) If the current flows in the direction of the angle @ from the magnetic 


south (Fig. 5), then the current density is expressed as 


J=cE (40) 
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where E is the component of E in the direc- 
tion of the current, and o the effective con- 
ductivity given by 


3 


Eastward 


Fi Fp xOyy tO py” 
 6,28iN*0+ ayycos?d 


(41) x 

For the N-S current (@=0° or 180°) and 
E-W current (9=90° or 270°), the effective 
conductivities become 


-_ Oxy? 7; 
o,=0 SF 2 
cet acl (42) i 
= o — Southward 
Be od 


respectively. 

The vertical distribution of effective conductivity is shown in Figs. 6A to 6D. 
From the Figures, it is noticed that the conductivity of the E-region is most predomi- 
nant in comparison with that of other regions. 

(iv) Following Martyn, we denote the height-integrated conductivity by 
>d\=Sedh). For convenience, we devide the ionosphere in three regions as follows: 

F-region; 400km to 170km 

E-region; 170km to 82km 

D-region; 82km to 60km. 
Assuming that in each region of these, there is no height gradient of horizontal 
electric field, then the total current in each region is given by 


T= (Jdh=§o,,dh-E,+§0,ydh-Ey=Di..E ct DieyEy (44) 
L,= (Jydh = —§a,ydh- E+ Soyydh- Ey=—\yE.t+ DyyEy (45) 
Let © be the angular distance to the direction of resultant current J having two 


components J, and J,, from the magnetic south, and E the component of E in the 
current direction, then, from (44) and (45) the following expression is obtained: 


I=XE (46) 
where 
Se 2 +2 Ps 2 +31, 2 
° o >3,,5in?0 133.com@ (47) 


which is the effective integrated conductivity. 
For the N-S and E-W directions, (47) becomes 


a +See (48) 
pi yy + (49) 


respectively. At the poles both these expressions simplify to 
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Figs. 6A to 6D—Altitude vs effective conductivity in different directions (6’s) at 
different magnetic latitudes (¢’s). 


a x2 

de= Dy = dh + -=2) (50) 
At the magnetic equator 

SL=s). Sy = o,dh=)};' (51) 


The numerical results of these integrated and effective integrated conductivities for 
above-mentioned three regions are shown in Tables 1 and 2, 
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Table 1—Values of height-integrated conductivities for different ionospheric regions 
expressed in e.m.u., and their ratios. 


Region 


= = = Zs SY | Sesh | Sw/as 
F 4.3010- | 8.8110-| 1.0710-"| 8.81x10-| 8.81x10- |] 0.012 1.00 
E 6.90 10-5 | 5.45x10-9 | 9.44x10-9 | 2.18x10-8 | 7.58x10-8 1.73 3.48 
D 7.78 X10-” | 3.58x10-! | 3.2610-| 6.54x10-%| 7.97x10- |] 0.91 1.22 


ee a 


Table 2—Values of height-integrated and effective height-integrated conductivities for 
different ionospheric regions at different magnetic latitudes expressed in e.m.u. 


Dip |Region| Sie De Sv 53H Sy 
F 4.30 10-3 0 8.81 10-1 || 4.30x10-3 |. 8.81x10- 
E 6.90 x 10-5 0 TERI] Ff 6.901027" 7.58 lO 
D 7.78 X10- 0 7297 SARE ET BOC 102 in 9 (eas? 
F 3.5210-9 | 2.14K10-%| 8.81x10- || 3.52x10-9 | 8.81 x 10-1 
E 2.16x10-% | 1.89x10-8 | 5.64x10-9 || 8.46x10-§ | 2.21x10-5 
D 5.47 X10- | 2.95x10-1 | 5.46x10-2 || 7.06x10- | 7.05 10-1” 
il L17 X10 | 22410" | S.81X10-2 17 X10 "8-81 X10 
E 72k X10°9 (CEOS K IOS [75.44 10 .2:90X10F 1) 22910) 
D -4.3.91X10-"| 2.99 10-© |-3,91 «x 10- |||. 6,19 x10-2| 6.19 x 10-2 
F 8.81x10-" | 1.07x10-#| 8.81x10-% | 8.81x10- | 8.81x10-% 
E 5.45x10-9 | 9.4410 | 5.45x10- || 2.1810-§ | 2.1810 
D 3.58 x 10-9 | 3.26 10-1? 6.54 x 10-10 


3.58 x 10-10 | 6.54 x 10-1 


6. Concluding remarks 

As mentioned above, the vertical distribution of electrical conductivity in the 
upper atmosphere depends on the atmospheric and.ionospheric models. Our calcula- 
tion is based on these models deduced mainly from the results of rocket measure- 
ments, therefore, this result may be no more than an example. But, it is asserted 
that even if the atmospheric or ionospheric models have a few modifications, the 
condvctivity of the E-region, especially near 100km, is most predominant, which is 
also pointed out by Hirono [14] and Martyn [15], so that this region is most probably 
effective to the electrical currents producing daily magnetic variations. This finding 
will be supported more strongly, if tidal motions in the F-region are reduced by 


induction drag as pointed out by Cowling [9, 16]. 
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LETTER TO THE EDITORS 


Remarks on Annual and Diurnal Variation of Cosmic Ray Intensity 


In previous paper (1), we estimated the coefficients of the temperature effect 
upon the cosmic ray data of Huancayo (12.°0S), Cheltenham (38.°7 N) and Godhavyn 
(69.°2 N) by comparing them with the corresponding aerological data such as San Juan 
(17° N), Washington (38.°9N) and Thule (77° N). And we found that the coefficients 
are larger for the data of lower latitude than higher latitude and confirmed that the 
amplitude of diurnal variation of cosmic ray intensity at Huancayo exceeds 0.6% or 
more if the data are corrected for the temperature effect. In these consideration, 
however, we had disregarded the positive temperature effect pointed by A. Duperier. 
Therefore, we have investigated the influence of this effect upon the annual and 
diurnal variations of cosmic ray intensity rather qualitatively. In what follows, we 
shall use the altitude difference between 100mb and 200 mb levels and the height of 
200 mb level above sea level as the expression of upper air temperature and the mean 
atmospheric temperature respectively. 


Altitude difference between 


4400 100 mb and 200 mb. Altitude difference between 


100 mb and 200 mb. 


Height of 200 mb. 


Height of 200 mb. 


116 122 
Jan. Feb. Mar. Apr. May, Jun. Jul. Aug. Sep. Oct. Nov. Dec. 02468 0 NU BBS Das 
Fig. 1 Annual variation of upper and lower Fig. 2 Diurnal variation of upper and lower 
atmospheric temperature (Tokyo). atmospheric temperature (Tokyo). 


Fig. 1 and Fig. 2 shows the annual and diurnal variations of upper and lower 
atmospheric temperature near Tokyo (35.°7N). The latter is averaged from bi-hourly 
radio-sonde data of 63 days which are obtained during 1944-1946 at Aerological 
Observatory Tateno, 60km north from Tokyo and the former is averaged from 
3 years (1949-51) radio-sonde data of Haneda Airport Station. These figures manifest 
that the annual variations of upper air temperature is opposite to that of lower 
atmosphere, namely, the upper air temperature decreases in summer owing to the 
increase of water vapours in lower atmosphere by which the heat radiation from the 
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earth’s surface is mitigated. On the other hand, in the case of the diurnal variation, 
the temperature of upper atmosphere increases with that of lower atmosphere by the 
insolation in day-time. Of course, the air is almost transparent for the solar radiation, 
hewever, we can expect that the higher the altitude, the larger the amplitude of 
diurnal variation of air temperature might be, because the thermal equilibrium in the 
upper atmosphere is wielded by the solar ultraviolet radiation and by its absorbers 
such as atomic Nitrogen, Oxygen and molecule of ozone abandant in the upper air (3). 
At present, the obtained upper air data are insufficient for the quantitative treatment 
of these problems. It is, however, quite obvious that the variations of cosmic ray 
intensity at the ground have the positive correlations with that of upper air tempera- 
ture in both cases of annual and diurnal variations. On the contrary, the mean 
atmospheric temperature and cosmic ray intensity are in the negative correlation in 
the case of the seasonal variation while they are positive in the case of the diurnal 
variation. Thus if we take the positive temperature effect into consideration, the 
values of the temperature coefficient due to the decay effect of mu-meson component 
becomes almost the same for each latitude as shown in Table I in which the tempe- 
rature coefficients are taken such values that the annual variations of cosmic ray 
intensity at sea level are cancelled. 


Table 1 
Amplitude of : sy 
ey difference pene . Went nee T Re sa 3 Temperature 
100 & “300 mb. pee sEvee Intensity. Coefficient. Coefficient. 
San Juan 80m 130m Huancayo 1.1% 0.008% /m —3.5% /km (—9.0% /km) 
Washington 150m 460m Cheltenham 2.6% * —3.1 {=5.590cee 
Thule (200-250 m) 1700m Godhayn 5.5% - —3.2~2.9 (—3.4 ,, ) 


The coefficients in last column indicate the values estimated without regarding 
to the positive temperature effect as in previous paper. The value of positive tempera- 
ture coefficient 0.008% /m is evaluated from A. Duperier’s result (2) 0.129¢/C by the 
following transformation 


aa. OL oh 
6h eT 
oh 
et i i 
then 7 afr 
273 
=0.122 "00 ==0-008 %/m 
because h=hy ( 1 +77) 


where /,=4200m is the average altitude differences between 100mb and 200 mb. levels 
for the middle latitude and T is mean temperature in °C. 


Before we make the’ similar considerations pertaining to the diurnal variation 
> 
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it would be better to note the magnitude of diurnal variation of cosmic ray intensity 
at various latitudes. Fig. 3 shows the diurnal variations of cosmic ray intensity 


Godhavn (69.°2N) 


Thule (77°N. summer) 


—_—— — . 
_ = — 
oat —_— _——— 


Tokyo (36°N) 


Ob 24, BUG 102 N14. Gets 20.225 24 
0 2 4 6 8 10 N 4 16 18 20 22 ma Fig. 4 Diurnal variation of upper air tem- 
perature indicated by the altitude 
difference between 100 mb and 200 mb 
levels at three different latitudes. 


Fig. 3 Diurnal variation of cosmic-ray inten- 
sity at different latitudes. 


obtained from the cosmic ray data from 1938 to 1946 at three different stations (4). 
From this figure, we can see that the amplitude of diurnal variation of cosmic ray 
decreases with latitude, though the statistical significance of these data are not 
testified. As for the aerological data, we find similar changes as shown in Fig. 4. 

Thus, applying the corrections of temperatureeffects upon the diurnal variation 
of cosmic ray, we get the results as shown in Table II. (As the aerological data for 
Cheltenham, we used that of Tokyo (35.7°N) owing to the lack of aerological data of 
Washington available to estimate the diurnal variation in upper atmosphere and this 
substitution may be allowable in our qualitative consideration.) 


Table I 
| Amplitude of Amplitude of cosmic ray intensity 
altitude differ- | Amplitude of 
ence between | 200 mb level. Corrected by yeaa parent dor phat 
| 100 & 200 mb. barometric effect. effect. perature effect. 
San Juan 60m 100m Huancayo 0.45% 0.80% 0.32% 
Tokyo 50m 60m Cheltenham 0.32% 0.51% 0.12% 
Thule (<20m) (<.30m) Godhavn 0.25% (0.35%) (0.1 %) 


The values in parenthesis mean the presumed one. Obviously, these estimations 
are not definitive from the statistical point of view but it seems likely that the 
influence of positive temperature effect upon the diurnal variation of cosmic ray 
intensity at the ground are not negligible and the amplitude of the diurnal variation 
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of cosmic ray intensity at sea level decreases with latitude even though the positive 
temperature effect is taken into account, as shown by dotted curves in Fig. 3. 
On these respects, it is necessary to correct much more precise data of cosmic 
ray and aerology. Furthermore, as regards to the coefficient of positive tem- 
perature effect, the cited value seems too large for the correction of the data of 
ionization chamber. To clarify these circumstances, we have tried the theoretical 
calculations by taking the energy spectrum into account and the results will be 
published in progress of theoretical physics in near future. 


(1) 
(2) 
(3) 
(4) 


(5) 


K. MAEDA 
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September 30. 
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